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Abstract: Neuronal activity serves as a gateway between external stimu- 
lus from environment and the brain, often inducing gene expression changes. 
Alternative splicing (AS) is a widespread mechanism of increasing the num- 
ber of transcripts produced from a single gene and has been shown to alter 
properties of neuronal genes, such as ion channels (Xie and Black 2001) 
and neurotransmitter receptors (Mu et al. 2003). Patterns of neural tissue- 
specific AS have been identified, often regulated by neuron-specific splicing 
factors that are essential for survival (Jensen et al. 2000; Li et al. 2007), 
demonstrating the importance of AS in neurons. In vitro studies of neu- 
ronal activity found AS changes in response to neuronal activity in addition 
to transcriptional ones, raising the question of whether such changes are 
recapitulated in vivo on behaviorally relevant timescales. We developed 
a paradigm for studying physiological neuronal activity through controlled 
stimulation of the olfactory bulb, and performed RNA-Seq transcriptome 
analysis of olfactory bulbs from odor-deprived and stimulated mice. We 
found that physiological stimulation induces large, rapid and reproducible 
changes in transcription in vivo, and that the activation of a core set of 
activity-regulated factors is recapitulated in an in vitro model of neuronal 
stimulation. However, physiological activity did not induce global changes 
in post-transcriptional mRNA processing, such as AS or alternative cleavage 
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and polyadenylation. In contrast, analysis of RNA-Seq from in vitro stim- 
ulation models showed rapid activity-dependent changes in both transcrip- 
tion and mRNA processing. Our results provide the first genome-wide look 
at neuronal activity-dependent mRNA processing and suggest that rapid 
changes in AS might not be the dominant form of transcriptome alterations 
that take place during olfactory rodent behavior. 



1 Introduction 

Neuronal activity links environmental signals and the brain, often induc- 
ing activity-dependent gene expression changes that are essential for synap- 
tic plasticity, wiring of the nervous system during development (Lin et al. 
2008), and learning and memory (Flavell and Greenberg 2008). Several 
activity-dependent protein-coding genes have been identified, such as the 
protein-synthesis independent Immediate Early Genes (IEGs) (Flavell and 
Greenberg 2008), and the molecular mechanisms controlling their expres- 
sion and activation of their targets have been studied intensively. Mouse 
knockout studies have shown that loss of IEGs like c-fos in specific brain re- 
gions results in neuronal cell death and interferes with learning and memory 
tasks, indicating that activity-dependent transcription is critical for brain 
function (Zhang et al. 2002). Furthermore, deletions and breakpoints in the 
activity-dependent transcriptional activator CREB-Binding Protein (CBP) 
cause Rubinstein-Taybi syndrome, underscoring the relevance of these fac- 
tors for human genetic disease (Petrij et al. 1995; Flavell and Greenberg 
2008). 

Post-transcriptional mechanisms of controlling gene expression, such al- 
ternative splicing (AS) and alternative cleavage and polyadenylation (APA) 
have also been shown to be regulated by neuronal activity (Li et al. 2007; 
Flavell et al. 2008). AS provides a way of diversifying the proteome and has 
been shown to alter properties of neuronal genes, such as ion channels and 
neurotransmitter receptors. Patterns of neural tissue-specific of AS have 
been described, often regulated by neuron-specific expression of splicing fac- 
tors (Li et. al., 2007). In vitro studies of neuronal activity documented 
post-transcriptional changes such as AS changes in response to neuronal ac- 
tivity in addition to transcriptional ones (Xie et. al. 2001, Ehlers et. al., Lee 
et. al. 2007) raising the question of whether such changes are recapitulated 
in vivo on behaviorally relevant timescales. 

Activity-dependent gene expression is typically studied in culture by ei- 



2 



Downloaded from http://biorxiv.org/on September 18, 2014 



ther KC1 stimulation (Lin et al. 2008; Kim et al. 2010), long-term incubation 
with ion channel antagonists like bicuculline and tetrodotoxins (Ehlers 2003; 
Mu et al. 2003) or neurotransmitter receptor agonists (e.g. NMDA recep- 
tor agonists (Li et al. 2007)), potentially inducing neuronal depolarization 
at levels that are not representative of physiological activity in an intact 
brain. While there have been genome-wide studies of activity-induced gene 
expression in vitro (Lin et al. 2008; Flavell et al. 2008), the prevalence of 
activity-dependent post-transcriptional changes in vivo is unknown. 

The olfactory bulb is an attractive model for studying neuronal activity- 
induced changes in vivo, since: (1) the synaptic pathway from olfactory 
sensory neurons in the nasal epithelium to the olfactory bulb is mapped 
and well-characterized, (2) the olfactory bulb is relatively homogeneous, 
easily dissectible brain region, made up of over 90% granule cells (Lledo 
et al. 2008), (3) olfaction is the primary sensory modality for rodents and 
is critical for survival, foraging and mating, and (4) the olfactory bulb is 
a highly conserved ancient brain region, similar in anatomy and function 
across vertebrates. 

Here we use an olfactory stimulation/deprivation paradigm and high- 
throughput mRNA sequencing (RNA-Seq) to measure genome- wide activity- 
dependent changes in the olfactory bulb transcriptome in vivo. We found 
that physiological stimulation induces large, rapid and reproducible changes 
in transcription, and that the induction of a core set of activity-regulated 
factors is recapitulated in an in vitro model of neuronal stimulation. How- 
ever, physiological neuronal activity did not result in global changes in post- 
transcriptional mRNA processing, such as AS or APA. In contrast, analysis 
of RNA-Seq from in vitro stimulation experiments shows rapid activity- 
dependent in both transcription and mRNA-processing. 

2 Results 

We developed a stimulation paradigm to induce neuronal activity in the ol- 
factory bulb using odorants (Figure 1A). Control mice were kept in a sterile 
cage with a constant supply of odorless air (pumped through an activated 
carbon filter, described in Methods) for 12 hrs and then immediately sacri- 
ficed. To induce neuronal activity in the olfactory bulb, mice were exposed 
to a cocktail of odorants (listed in Supplementary Table 1) chosen to max- 
imize activity across the olfactory bulb based on their empirical glumeruli 
activation profiles (Johnson and Leon 2007). 

Two stimulation time courses were used to induce neuronal activity. The 
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first wave of transcriptome changes was measured in a short time course, 
where animals were placed in the odor ant- free cage for 12 hrs and then 
received 30 min of odorant stimulation. A second time course was used 
to measure effects downstream of the first wave and changes with slower 
kinetics, where mice received the same treatment as the 30 min animals but 
were given an additional 90 minutes of clean air ("120 min" time course.) 
(Figure IB). Mice in the stimulated conditions exhibited highly exploratory 
behavior and movement upon induction of the odorants and throughout its 
duration, suggesting that animals were not desensitized to odorants during 
the experiment. 

We prepared RNA-Seq libraries from whole olfactory bulbs (see Meth- 
ods), obtaining a total of ~320 M uniquely mapping reads from 9 mice. 
Reads were mapped to the genome and splice junctions (see Methods) and 
then used to compute gene expression profiles for control and stimulated 
mice. Immediate early genes (IEGs), such as c-fos, were rapidly induced in 
all 30 min stimulated mice, with upto ~22 fold change in expression rela- 
tive to controls, shown in (Figure 1C). Induction of the IEG Egrl mRNA 
was confirmed by in situ hybridization (Supplementary Figure SI), corrob- 
orating RNA-Seq results, while induction of c-fos protein was confirmed by 
immunohistochemistry one hour following stimulation (data not shown.) 

Hierarchical clustering of gene expression profiles for all mice revealed 
that replicate samples cluster together and exhibit high correlations in gene 
expression (Figure ID). Furthermore, 120 min stimulated mice were more 
similar globally to the control mice than to 30 min stimulated mice, reflecting 
the fact that most activity-induced gene expression changes are transient and 
return to near-baseline levels by 120 min (Figure ID). 

To determine the overall kinetics of transcriptional changes, we per- 
formed hierarchical clustering analysis of gene expression profiles across the 
control and stimulation time points (Figure 2A). Several clusters of genes 
emerged that vary in the direction and shape of expression profiles. IEGs 
such as c-fos and Egrl fit the expected "fast rise, fast fall" genes cluster, 
where expression is quickly upregulated but returns to baseline levels by 
120 min (Figure 2A, Cluster 8). A distinct cluster of genes showing stable 
upregulation was also detected, and contained IEGs with slower kinetics like 
Homer 1 (Figure 2A, Cluster 5). 

Clusters were classified into three general categories: clusters where the 
net effect of activity is to upregulate expression ("Enhanced" genes, whose 
expression is increased in both stimulated conditions relative to controls), 
genes where the net effect of activity is to downregulate expression ("Re- 
pressed" genes, whose expression is decreased in both stimulated conditions), 
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and genes where the effect of activity switches between the two stimulation 
conditions relative to control ( "Biphasic" genes) . This demonstrates hetero- 
geneity in the effect of activity on transcription. 

Known IEGs were induced reproducibly in all stimulated mice compared 
to controls (Figure 2B) and exhibited the expected kinetics. We also ob- 
served increased transcription of several activity-dependent factors impli- 
cated in neurogenesis/cell-cycle or enriched in newborn cells (Supplemen- 
tary Figure S2), such as Reelin (Hack et al. 2002) and TenascinR (the lat- 
ter known to be activity-dependent (Saghatelyan et al. 2004)), suggesting 
that whole olfactory bulb RNA-Seq provides sufficient sensitivity to detect 
changes in genes expressed in specific cell subtypes. 

Recently, high-resolution RNA-Seq and ChlP-Seq was performed on 
mouse neurons stimulated in vitro look at activity-regulated RNAs and 
transcription factors controlling their expression (Kim et al. 2010). Cul- 
tured neurons were profiled following 1 and 6 hours of KC1 treatment. 

We applied the computational methodology used to analyze olfactory 
bulb RNA-Seq data to examine activity-regulated transcription and mRNA- 
processing in KCl-treated and control cultured neurons. Fold changes in 
expression between the first stimulated conditions and controls in the olfac- 
tory bulb and neuronal cultures were compared (using 30 min stimulation 
and OB control versus 1 hr KC1 stimulation and KC1 control, respectively, 
Figure 2C). As expected, the majority of genes show fold changes < 1.5-fold 
upon stimulation in both experiments. Strikingly, a group of genes were 
strongly upregulated — with fold changes upwards of 2-fold — in the stimu- 
lation conditions of both experiments (red points, Figure 2C). This group 
comprised the genes with the largest fold changes overall, though the mag- 
nitude of fold changes was consistently larger in neuronal culture. ChlP-Seq 
analysis showed that the majority of these commonly induced factors were 
bound by Npas4, consistent with the possibility that the maintenance of 
their expression is regulated in part by Npas4 (data not shown.) We also 
noted a smaller number of genes that are downregulated with activity in 
both experiments. 

Several novel predicted protein-coding genes were upregulated in both 
experiments (data not shown), while the noncoding small nuclear RNA 7SK, 
implicated in transcriptional control (Peterlin et al. 2012), was downregu- 
lated in both experiments. 7SK is traditionally thought to be transcribed 
by RNA Pol III, though recent evidence suggests it associates with RNA 
Pol II and has a histone signature of Pol Il-transcribed loci (Barski, 2010). 
This suggests that regulatory noncoding RNAs also play a role in activity- 
dependent regulation in neurons. 
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We next asked whether activity induces changes in mRNA-processing. 
Several activity-dependent skipped exon events and their associated cis- 
regulatory elements have been identified (Li et al. 2007; Xie et al. 2005). 
It was also previously shown using exon arrays that long-term visual de- 
privation of mice followed by light stimulation induces changes in alterna- 
tive cleavage and polyadenylation of a subset of transcripts in visual cor- 
tex (Flavell et al. 2008). Therefore, we focused on alternatively skipped 
exons and alternative cleavage and polyadenylation sites as representative 
instances of mRNA-processing events. 

We found that activity did not induce global changes in olfactory bulb 
splicing patterns, and similar trends were observed for alternative cleavage 
and polyadenylation events (data not shown.) By contrast, the same analy- 
sis revealed strong and rapid activity-dependent splicing changes in neuronal 
cultures. The KCl-regulated exons we discovered corroborated findings of 
activity-dependent splicing in culture (Li et al. 2007; Xie et al. 2005). Cer- 
tain exons showed strong and rapid switch-like changes in splicing within 
1 hr of KC1 treatment (Figure 3A,B), which is significantly quicker than 
the 24-48 hr incubations with depolarizing agents used for previous in vitro 
measurements of activity-dependent splicing (Lee et. al., 2007). Some exons 
were also strongly regulated by activity in the 6 hr KC1 time point. The 
Adyacaplrl gene contains a calcium responsive RNA element (CaRRE) in 
the intron upstream of Exon 14 (Figure 3B) that causes repression of the 
exon upon depolarization in CaMK IV-dependent manner (Li et al. 2007). 
We found that this exon is strongly skipped with activity and is near absent 
in the 6 hr KC1 condition (Figure 3C). Similarly, Arghgap21 contains an 
exon that decreases ~3-fold in inclusion levels by 6 hours of KC1 treatment 
(Supplementary Figure S2). 

3 Discussion 

Our study provides an in vivo genome-wide look into physiological activity- 
dependent mRNA processing in the mouse olfactory bulb transcriptome. 
Our results suggest that rapid transcriptional regulation, but not mRNA- 
processing, predominates the olfactory bulb transcriptome response to neu- 
ronal activity. This has implications for the study of gene expression changes 
induced by olfaction-based behavioral experiments in rodents. 

A systematic comparison of our in vivo stimulation of olfactory bulbs 
and a commonly used in vitro model of neuronal stimulation suggests that 
a core set of activity-regulated factors is shared across diverse neuronal cell 
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types, but that mRNA-processing changes are not. Recent work shows that 
not all brain regions employ a high degree of alternative splicing. For exam- 
ple, recent study of cerebellum (Pal et. al., 2011) showed that throughout 
development of the brain region, alternative promoter usage is a much more 
commonly used form of generating gene product diversity, compared with 
alternative splicing. Both the cerebellum and the olfactory bulb are ancient 
and highly conserved brain regions in vertebrates, leaving open the possi- 
bility that evolutionarily newer regions (like cortex) do employ alternative 
splicing and mRNA-processing more widely. 

We cannot exclude the possibility that some of the discordance ob- 
served between our in vivo model and in vitro stimulation models are due 
to the difference in developmental stages of the cells (neuronal cultures 
originating from embryo, our tissues from an adult animal.) In principle, 
there could be developmental stage-specific activity-dependent instances of 
mRNA-processing. Future studies using transcriptome sequencing in of ac- 
tivated neurons in distinct developmental stages and brain regions could 
resolve these issues. 

4 Methods 

4.1 Olfactory stimulation and deprivation setup 

C57BL6 mice (Charles River laboratory) ages 8-11 weeks were placed in 
an autoclaved ethanol-sterilized cage and placed in a clean fume hood for 
12 hours overnight, with a steady supply air filtered by a compressed air 
activated-carbon filter (Hankison, Inc.). Filter was hooked into laboratory 
compressed air nozzle mediated by an air pressure gauge. A timer-controlled 
light source programmed to the animal housing facility's light/dark cycle 
was placed above the cage during experiments. All animals were sacrificed 
between 10 am and 12 pm. 

4.2 Odorants preparation 

A set of 12 odorants (listed in Supplementary Table 1) was selected based 
on each odor ant's activation profile of specific glomeruli or subregions of the 
olfactory bulb as described in (Johnson and Leon 2007) and related stud- 
ies. Odorant mixtures were comprised both glomerulus-specific and broad 
activators of the olfactory bulb (such as valeric acid and methyl valerate), 
with many of the odorants having food- associated scents. Odorants were 
diluted in Light Mineral Oil (Sigma) and pipetted onto damp pieces of filter 
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inside a capped Erlenmeyer flask, which was plugged into the air supply of 
the mouse cage during stimulation and pumped in at 3 psi. 

4.3 Tissue preparation and extraction of RNA 

From each animal, one hemisphere's olfactory bulb (OB) was used for RNA 
extraction and RNA-Seq library preparation and the remaining bulb taken 
for RNA in situ hybridization and immunohistochemistry to confirm in- 
duction of immediate early genes Egrl and c-fos. OBs for RNA-Seq were 
dissociated in Trizol (Invitrogen) and RNA extracted using standard Trizol 
extraction, followed by a run on a BioAnalyzer to profile RNA quality and 
obtain RNA Integrity Number (RIN). The remaining intact OBs (attached 
to forebrain) were prepared for sectioning by immersion in TissueTek in- 
side a plastic well, and then flash frozen in crushed dry ice and isopropyl 
alcohol. Cryostat coronal sections of frozen OBs were prepared for in situ 
hybridization. 

4.4 Sample preparation for RNA-Seq 

RNA-Seq libraries were prepared based on protocols described elsewhere 
(Mortazavi et al. 2008). Briefly, two rounds of polyA+ RNA selection were 
performed using poly(T) beads. Remaining polyA+ RNA was fragmented 
and prepared for Illumina sequencing: paired-end adaptors were ligated onto 
cDNAs primed with random hexamers and products were size-selected by 
agarose gel electrophoresis. Bands ranging from ~250-300 bp were cut out 
using GeneCatcher gel excision tips and purified. Samples were PCR am- 
plified for 13 cycles prior to sequencing. 

4.5 Mapping and computational analysis of RNA-Seq data 

Reads were mapped to mouse genome (mm9) and a set of splice junctions 
(Wang et al. 2008) using bowtie (Langmead et al. 2009). Only uniquely 
mapped reads were used for downstream analyses. Quantitation of exon 
abundances was performed using MISO (Katz et al. 2010). Gene expression 
levels quantitated by computing RPKM values (Mortazavi et al. 2008) using 
only reads from constitutive exons. 

SOLiD mRNA-Seq sequencing data from (Kim et. al., 2010) was ob- 
tained from GEO (Accession code: GSE21161) and mapped to the same 
mouse (mm9) indexed genome plus splice junctions, built using the bowtie 
color index (-C) option. 



8 



Downloaded from http://biorxiv.org/on September 18, 2014 



Expression profiles of genes (in RPKM units) showing 1.5-fold change 
or greater between any pair of olfactory bulbs were used to perform hier- 
archical clustering using the Euclidean distance metric in Figure ID. These 
expression profiles were z-score normalized to perform hierarchical cluster- 
ing analysis of each gene's expression profile over the three olfactory bulb 
conditions (Figure 2A). All genes expression analyses and clustering were 
performed using custom-built scripts written in the Python programming 
language, using scipy/numpy packages. 

4.6 Data availability 

mRNA-Seq data from this study are available at GEO (Accession code: 
GSE58843). 
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Figure 1: Olfactory stimulation paradigm followed by RNA-Seq. 
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List of Figures 

1 (A) Adult (8-11 weeks old) mice were kept in odor- free en- 
vironment for 12 hrs using an activated-carbon air filter and 
then stimulated with a mixture of odorants (odor-deprived 
animals used as controls.) Immediately following stimulation 
or deprivation mice were sacrificed, and olfactory bulbs ex- 
tracted from which RNA-Seq libraries were prepared. (B) 
Timelines of olfactory deprivation/stimulation experiments. 
Control animals were odor-deprived for 12 hours. Stimulated 
animals were odor-deprived for 12 hours, followed by either 
exposure to 30 minutes of odorant stimulation only ( "30 min" 
time course) or 30 minutes of stimulation followed by 90 min- 
utes of clean air ("120 min" time course). (C) RNA-Seq reads 
mapping to the c-fos immediate early gene from 9 mice (3 
mice per condition) used in the experiment. Control condi- 
tion replicates shown in blue, 30 min. condition in red and 120 
min condition. Read densities correspond to number of reads 
overlapping a position (log 10 units) and arcs correspond to 
exon-exon junctions whose width is proportional to number 
of junction reads. (D) Average-linkage hierarchical cluster- 
ing on genes with fold change in expression greater than X 
between any pair of time points. Replicate mice cluster to- 
gether, while the control gene expression profiles are globally 
more similar to 120 min mice than 30 min mice 10 
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2 (A) Hierarchical clustering of gene expression profiles (av- 
erage of replicates) across time points. Mean z-score nor- 
malized expression values are plotted along with shaded re- 
gions around standard deviation of scores. Mean expression 
of genes whose expression is enhanced by activity in both 
time points is shown in green ("Enhanced" genes), those 
whose expression is repressed in red ( "Repressed" genes) , and 
genes that change direction between the two time points in 
response to activity (relative to the control) in black ( "Bipha- 
sic" genes.) Numbers in parentheses on y-axis of each cluster 
indicate the number of genes in that cluster. (B) Gene ex- 
pression changes from olfactory bulb RNA-Seq for immediate 
early genes c-fos, Arc, Egrl, Npas4 and Homerl. Expression 
levels measured in 'reads per kilobase per million mapped 
reads' (RPKM) in controls, 30 min stimulation and 120 min 
stimulation are plotted in blue, red and orange, respectively. 
(C) Scatter plot of fold changes (in log 2 ) in expression be- 
tween OB 30 min stimulation time point and OB control ver- 
sus neuronal culture 1 hr KC1 and culture controls. Grey 
points (plotted with alpha transparency) correspond to genes 
with fold changes smaller than 1.5 in either experiment, red to 
genes upregulated with fold changes > 1.5-fold in both exper- 
iments, and black to genes downregulated with fold changes 
> 1.5-fold in both experiments. Select commonly upregulated 
IEGs and downregulated 7SK noncoding RNA are labeled. . 11 

3 Effect of activity on mRNA-processing. (A) Rapid KCl-induced 
repression of alternative exon in Arpp21 gene. Read densities 
(in base log 10 ) along exons are plotted along with junctions 
as arcs. By 1 hr KC1, the alternative exon (marked by arrow) 
is nearly absent from transcripts. (B) Exon (marked by ar- 
row) in Adcyaplrl is strongly skipped in 6 hr KC1 condition, 
but not 1 hr KC1, relative to control. Upstream intron of 
repressed contains CaRRE (sequence shown in capital bold 
letters, marked by curved arrow) known to inhibit splicing in 
depolarization-dependent manner 12 

SI (A) mRNA fluorescent in situ hybridization for Egrl in coro- 
nal olfactory bulb sections from control ("clean air") and 30 
min. stimulated mice. (B) RNA-Seq estimates of gene ex- 
pression (in RPKM) for Egrl in controls and stimulated mice. 13 
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RNA-Seq expression levels for genes either associated with 
adult neurogenesis or enriched in newborn/migrating cells 
(Reelin, TenascinR, p21/Cdknla, Btg2) that show activity- 
dependent regulation in olfactory bulb samples 
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